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Abstract. We construct singular solutions to special Lagrangian equa- 
tions with subcritical phases and minimal surface systems. A priori 
estimate breaking families of smooth solutions are also produced cor- 
respondingly. A priori estimates for special Lagrangian equations with 
certain convexity are largely known by now. 



1. Introduction 

In this paper, we construct singular solutions to the special Lagrangian 
equation 

n 

(1.1) arctan = 

1=1 

with subcritical phase |0| < {n — 2) tt/2, where Xt are the eigenvalues of 
D'^u, and the minimal surface system for fc-vector valued functions of n- 
variables 

n ^ 

(1-2) ^gU=Y, —a., iV99''d,^U) = 0, 



where the induced metric 



I+{DU f DU. 



Equation (jl.ip is the potential equation for (jl.2p with solutions U = Du. 
The Lagrangian graph {x,Du{x)) C M" x M" is called special and in fact 
volume minimizing when the phase or the argument of the complex num- 
ber (l -|- \/—lXi) • • • (l -|- \/—lXn) is constant 6, or equivalently u satisfies 
equation (jl.ip : see the work [HLl, Theorem 2.3, Proposition 2.17] by Har- 
vey and Lawson. The phase (n — 2)7r/2 is said critical because the level 
set { A G M"| A satisfying (jl.ip j is convex only when |0| > (n — 2) 7r/2 [Y2, 
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Lemma 2.1]. In dimension three, when |0| = 7r/2 or |0| = 0, vr, equation 
(|l.ip also takes the quadratic and cubic algebraic forms respectively 

(1.3) cj2(Z)2n) =AiA2 + A2A3 + A3Ai = l 
or 

(1.4) Au = del D'^u. 
We state our first main result. 

Theorem 1.1. There exist C^'^/C^'^^^) Cm = 2,3,4, ••• j viscosity solu- 
tions to with n = 3 and each Q G (~f > f ) i such that n™ E 
(ji,i/(2m-i) p (joo (Bi\{0}) for Bi C hut ^ C^'-^ for any 6 > 
l/(2m - 1) . 

Rotating forth and back, we obtain our second ("smooth") result. 

Theorem 1.2. There exist a family of smooth solutions to in Bi C 
with n = 3 and each fixed Q £ (~f ' f ) such that 

- but I Z:'^n^ (0) I oo as e ^ 0. 

For each it^ with small e fixed in Theorem 1.2, the Hessian |D^u^(0)| 
(in the max eigenvalue norm) is strictly larger that its nearby values in the 
three dimensional domain of the solution to a now uniformly elliptic equa- 
tion (II. ip . (It can be seen by Property 2.4 in Section 2 and tracing the 
eigenvalue dependency in Section 4.) This violates the maximum princi- 
ple. In contrast to the two dimensional fully nonlinear uniformly elliptic 
equations, it is classically known that the Hessian of any solution enjoys the 
maximum principle. To the solutions in the above two theorems, by adding 
quadratics of extra variables in higher dimensions n > 4, we immediately 
get the corresponding counterexamples for (jl.ip with all subcritical phases 
|0| < (n — 2) 7r/2. Furthermore, we convert our counterexamples to the ones 
for minimal surface system (jl.2p . 

Theorem 1.3. There exist a family of weak solutions U'"^ to in Bi C 

with n = 3, k = 3, and m = 2, 3, 4, • • • such that 

jjm ^ ^i,p X p ^ 2m + I ^„ mi) . 

2m — 2 

Furthermore, there exist a family of smooth solutions to il.S^) in Bi C M'^ 
with n = 3 and k = 3 such that 

\\U'\\lo.(^b^) <C but \DU' {0)\^oo ase^O. 

The vector valued functions C/™ are taken as Du"^ with u™ from The- 
orem 1.1, thus the first part of the theorem gives a negative answer to 
Nadirashvili's question whether there is an e improvement of W'^'^ solutions 
to special Lagrangian equation (II. ip in general. We are grateful for this 
question. In terms of minimal surface system (jl.2p . the question would be 
whether there is an e improvement of W^'^ solutions. 
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For special Lagrangian equation (ll.ip with critical and supercritical phases 
I ©I ^ (n — 2)7r/2 in dimension two and three, with very large phase |0| > 
(n — 1) 7r/2 in general dimensions, a priori Hessian and gradient estimates, 
and consequently, armed with the solvability of the Dirichlet problem with 
smooth boundary data to the now convex special Lagrangian equation (jl.ip 
in the critical and supercritical phase cases, the regularity of viscosity 
solutions were derived in [WYl] [WY2] [WY3] [CWY]. In passing, we also 
mention that the existence and uniqueness of the viscosity solution for the 
Dirichlet problem to strictly elliptic equation (II. 1[) is known (cf. [CWY, p. 
594]). In recent years, there has been a new approach toward the existence 
and uniqueness of viscosity or weak solution for the Dirichlet problem 
to strictly elliptic as well as degenerate elliptic fully nonlinear equations by 
Harvey and Lawson [HL2] [HL3]. 

Recently Nadirashvili and Viaduct [NV] constructed beautiful C^'^^^ sin- 
gular viscosity solutions to (jl.ip with subcritical phases |0| < tt/2 in di- 
mension three, relying on "brutal force" calculations (for the approximate 
solutions) and a hard and deep topological result in [EL] (for the injectivity 
of the gradient maps). 

For minimal surface equations, namely (jl.2p with k = 1, the gradient 
estimate in terms of the height of the minimal surfaces, is the classic result 
by Bombieri-De Giorgi-Miranda [BDM] , from which it follows the regularity 
of weak or viscosity solutions. For smooth solutions to (|1.2p with n = 2, 
Gregori [G] extended Heinz's Jacobian estimate to get a gradient bound in 
terms of the heights of the two dimensional minimal surfaces with any codi- 
mension. For smooth solutions to general minimal surface system (ll.2p with 
certain constraints on the gradients themselves, a gradient estimate was ob- 
tained by Wang [W], using an integral method developed for codimension 
one minimal graphs. Nonetheless, there do exist singular W^''^~ weak so- 
lutions to (]1.2p with n = 2; see Osserman [O]. Now gradient estimates for 
(jl.2p with k = 2 and n > 3 still remain mysterious and challenging. 

Our construction goes as follows. In the first stage, we solve the special 
Lagrangian equation (jl.ip with the critical phase by Cauchy-Kowalevskaya. 
The approximate solutions or initial data for the relatively "easier" corre- 
sponding quadratic equation (jl.Sp are built up via a systematic procedure, 
which allows us to have the approximation at arbitrarily high order (Prop- 
erty 2.1 and 2.2), and eventually those highly ("oddly" C^'^/('^"^~^)) singular 
solutions in Theorem 1.1 and Theorem 1.3. In the second stage, we take an 
"inversion" ^ rotation of the solutions from the first stage to obtain those 
singular solutions with phase (Proposition 3.1). The singular solutions 
with other subcritical phases are achieved via a preliminary "horizontal" 
rotation before the "inversion" ^ rotation (Step 1 of Section 3). Some re- 
marks are in order. Those U (n) rotations are "obvious" to produce for the 
U (n) invariant special Lagrangian equation (jl.ip . But it is by no means 
easy to justify that the special Lagrangian submanifold is still a graph in 
the rotated new coordinate system, thus a valid equation (jl.ip to work on. 
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(Earlier development of those U (n) rotations for (jl.ip can be found in [Y2] 
[Y3].) Here our elementary analytic justification for the "inversion" ^ rota- 
tion (Proposition 3.1) avoids a hard and deep topological formula of [EL], 
which was employed in [NV] . Lastly we point out that the Legendre transfor- 
mation (usually used for convex functions), is just the "inversion" ^ rotation 
followed by a conjugation for converting "gradient" graph (x, Du (x)) to the 
one [Du* [y) ,y) (now with saddle potentials u and u*). In the third stage, 
we kick in a little bit extra to the preliminary "horizontal" rotations of Stage 
2, then after the same "inversion" ^ rotation, we make up a corresponding 
little bit "backward" rotation to finally generate the desired family of smooth 
solutions in Theorem 1.2, which break a priori Hessian estimates for special 
Lagrangian equation (jl.ip with subcritical phase. Note that here one cannot 
produce those a priori estimate breaking family of smooth solutions by the 
usual way, that is to solve the Dirichlet problem with smooth approximate 
boundary data of the boundary value of a singular solution, as Theorem 
1.1 shows the non-solvability of smooth solution to the Dirichlet problem to 
(jl.ip of subcritical phase even with smooth boundary data. The Dirichlet 
problem to the saddle branch of (jl.4p or the equivalent (jl.ip with n = 3 and 
= was "invited" by Caffarelli, Nirenberg, and Spruck in [CNS]. 

In closing, we point out that any further regularity beyond continuity for 
continuous viscosity solutions to general special Lagrangian equation (jl.ip 
is unknown. We are also curious to know whether there exist other C^'" (no 
better) singular solutions to (jl.ip with, in particular, irrational exponents 
a between those odd reciprocals 1/ (2m — 1) . Meanwhile, we guess that all 
C^'" for a > I solutions to special Lagrangian equation (jl.ip with n = 3 
should be regular (analytic). This regularity for C^'^ solutions to p.ip in 
dimension three was shown in [Yl]. Earlier on, Urbas [U, Theorem 1.1] 
proved the regularity for better than Pogorelov solutions, namely all C^'" 
for a > 1 — ^ (convex) solutions to the (dual) Monge-Ampere equation 
In det D^u = In Ai -|- • • • -|- In = c are C^'^ and eventually analytic. 



2. CAUCHY-KOWALEVSKAYA with CRITICAL PHASE = f 

As a preparation for the constructions in the next three sections, we solve 
the following special Lagrangian equation with critical phase in dimension 
three by Cauchy-Kowalevskaya. The quadratic nature of the equation at the 
critical phase is easier to work with than the cubic nature of the equations 
otherwise. 

Our approximate solution P (x) to the equation 



(2.1) 



CT2 {D^U) = 

U3 (xi,X2,0) 
U (xi,X2,0) : 



{Auf - \D^u\ 
Ps (a;i,X2,0) 

P (xi,X2,0) 



1 or 



arctan Aj = ^ 
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is a polynomial of degree 2m 



1 2 ''^ 



coefficients e and ajs arc to be dc- 



where Z = xi + v — 1x2 = pexp [^y—U 
termined later. We construct this P satisfying the following four properties, 
so does u then, for |x| = r < Vm with positive depending only on m. 
Property 2.1. (T2 (^D^P^ — 1 = , here [r^] represents an analytic 

function starting from order k. Then the solution u coincide with P up to 
order 3m — 2 (> 2m for m > 2, 3, 4, • • • ). 

Property 2.2. The three eigenvalues of D^P, then also D^u satisfy 



A 



l+[ 



m— 11 



„2m— 2 ^ \ ^ r I ^\ J2m—2 



A2 = 1 + [r- 



-52 (m) r^™-^ < As < -6i (m) r' 
Property 2.3. The "gradient" graph 

(x, Du) = 

x,xi+0 ip) [r™-i] + [r^™] , X2 + [p) [r™"!] + [r^™] , 
ReZ™ + ^p2m-2^^ _ 2mexf^-^ + ep^ [r^^-S] + [^2™] 

Property 2.4. The gradient Du satisfies 

h {jn) r^""-^ < \Du {x)\ < 64 (m) r. 

We first find the equation near a quadratic solution. Let 



u = - {fJ,ixl + p2xl + nsxl) +w{s 



Then 



a2{D\)-l = ^[iAuf-\D' 



- 1 



1=1 



- 1 



= m {Aw - wu) + p,2 {Aw - W22) + P3 (Aw - Was) + ^ 

+P1P2 + li2P-Z + - 1- 

Set /xi = /X2 = 1 and /X3 = 0, we get 

2 



fA-u;)^ - \D^w\ 



02 {D^u) — 1 = Wii + U;22 + 2^33 + - {Aw} 



ID^wl 



Aw + l 



{Awf 



\D^w\ 
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where A = du + 822 + 2933. To make the right hand side of the above 
equation vanish at high orders, we choose w = h + Q + H, where 

h = Re Z™' 2:3, an ad hoc "harmonic" function; 
,2 



Q 



H = e 



-P 



^2m 
"^3 



x§, to match CJ2 {D /i) ; 

, to make eigenvalue A3 negative. 



Then 



CJ2 {D^u) - 1 = Ah +AQ + AH + ^ [{^hf - \D^hf + [r^^-Sj 



A simple calculation leads to 
Re [m (m - 1) Z""-^] 

D'h 



X3 



Im 
Re 



m [m - 1) 



X3 Re mZ"^ ^ 
X3 —ImmZ"^^^ 




It follows that 



C72 {D^h) = - [m (m - 1) p™-^] - m~ p 



2 „2m-2 



Thus 



XI + n?p''^-^ + CJ2 {D^h) = 0. 



AQ + a2 {D^h) = [m{m- 1) p"^-^] ^ xj 
Finally we fix the "harmonic" H satisfying AH = with 



ao 



2 • (2m - 2j + 2) (2m - 2j + 1) 



-a,'_i 



_ .+1 2i2m(2m-l)---(2m-2j + l) . 
^ ' 2242... (2jf IorJ-i> 

and e is still pending. Therefore, P = \ {x\ + x"^ + h + Q + H, satisfies 

£72 (D^P) - 1 = [r3™-3] . 

Now the analytic solution u to (j2.ip with initial data P follows from 
Cauchy-Kowalevskaya. As in [NV], considering the linear equation for dif- 
ference u — P, the Cauchy-Kowalevskaya procedure implies that the solution 
u coincides with P up to order 3m — 2 (> 2m for m > 2). Thus Property 
2.1 is verified. 

We move to Property 2.2. We have 



(2.2) D-^u 



1+ [1 



ImrnZ™-! + [1 

m'' .2m— 2 i u \ \r, 
— P + -"33 + 



,2m-21 



,2m-ll 
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Because the eigenvalues are Lipschitz functions of the matrix entries, we get 

Ai = 1 + [r^-^] 
A2 = 1 + [r^-^] . 

By the quadratic Taylor expansion of the isolated eigenvalue A3 in terms of 
the matrix entries near D^u (0) , we obtain 



A3 = «33 - uls - vm + 



p2m—2 _j_ |'j,2m-l 



for m > 2 

2m-ll 



c- r Orr, (Orr, 1 ^ ^2m— 2 , ~ 2m— 4: „2 I _|_ ~ ^2m— 21 m? Jim— 2 

e \—zmyzm — L) +022:3 p + ■ ■ ■ + am-ip J — yP 



+ [r 



H: 



33 



^ ^2m-2 

2 

The "harmonic" function H-^^ cannot have a definite sign near the origin, 
but with the help of — ^p^"*"^ and small e, we make A3 negative. Let rj be 
a small positive constant to be chosen shortly. 
Case 1: ri\x3\ > p. We have 



[-2m (2m - 1) 



.2m-2 ^ a2xl"'-'^p'^ + •••]=- [2m (2m - 1) + O (1) r/^ 



Note r/-\/l + 7/2 < \xs\ < r, then 



£ [2m (2m - 1) + O (1) ry^] | ^^^-2 < A3 < -e 



"I" 2 



2m(2)»-l)-0(l)i;^ 

/ ; \ 2m — 2 



„2m-2 



„2m-2 



Case 2: 77 |a;3| < p. Note rri/\/l + r/^ < p < r, we have 



[-2m (2m - 1) xf"-^ + aaxf^^V^ + 



,2m-2 " ' 



then 



and 



A3 



r 2 
m 



™^2m-2 



^2m — 2 

+o{l) 



„2m-2 



< A3 < <^ 



£0(1) 



> r 



2m-2 



+o(l) 

We first choose rj = r] (m) > small, next e = e (ry, m) > smaller, then 
there exist Si = 5 {rj, m) > and S2 = 62 (m) > such that 



for r < Vm- Here is within the valid radius for the Cauchy-Kowalevskaya 
solution u. 

Property 2.3 follows from u = P + [r^"*~^] . 
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Finally we prove Property 2.4. The upper bound is straightforward. For 
the lower bound, from Property 2.3, we have 



„mi\2\2 



m 



Case 1: x| > p. From = + a;| < (x| + l) we know 

la^al > r. 

Note that the other terms than —2mexf^~^ in (x) have the following 
asymptotic behavior near the origin 



x 



2m 
3 ' 



^ 2m — 2 



^ I |4m-3 



0{x] 



,2m \ 
3 J 



It follows that 



\Du{x)f > \u3{x)f = [-2mexf^-^ + 0{xf^)]^ 
> ,53 (m) > ,53 (m) r2(2— 1) 

for |x| < Vm with positive rm and ,53 (m) to be fixed shortly. 
Case 2: x| < p. From = p^ + x| < (p + 1) p, we know 



p > 



Then 



\Du{x)\'^ > u\ (x) + ul (x) = p2 + 2a;i[r'"] + 2x2[r"^] + 2[r 



mi2 



= p2 + o(p— 

- 2^^ - 2 - 



2(2m-l) 



for p < r < with the positive Vm to be fixed next. 

Now we choose positive ,^3 (rn) small and the small positive within the 
valid radius for Cauchy-Kowalcvskaya solution u and Property 2.2, Property 
2.4 is then completely justified. 

Since u (rmx) /r^ is still a solution to a2 [D'^u) = 1 in Si C i?^. We may 
assume the above constructed solution is already defined in Bi C . Note 
that D [u {r„ix) /r^] = Du (r^x) /r^ and [u (r„,,T) /r^] = D'^u (rmx) , 
we see that Property 2.2 and Property 2.4 are still valid in Bi with 5i, 62, S3 
replaced by r^~^Si, r^~'^62, r^~^62 respectively, and ,54 unchanged. 
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3. Rotate to subcritical phases |6| < f : proof of Theorem 1.1 

In this section, we carry out the construction of the singular solutions in 
Theorem 1.1 by "horizontally" and 7r/2 rotating the Cauchy-Kowalevskaya 
solutions from Section 2. The latter rotation, Proposition 3.1 is pivotal. 

Step 1. Let a £ [0, 7r/4). We will take a = 6/2 for 6 G [0, 7r/2) in Step 3 of 
this section. We make a U (3) rotation in : z' = e'^^'^z' and z^ = z^ with 
z = {z',zs) = {x',X3) + ^/^{y',y3) and z = {z',zs) = {x' ,xs)+^/^ {y' ^y^) . 
Because U (3) rotations preserve the length and complex structure, 9JT = 
{x,Dv{x)) for X £ Bi is still a special Lagrangian submanifold in the new 
coordinate system with parameterization 



(3.1) 



X = (xi cos a + ui (x) sin a, X2 cos a + U2 (x) sin a, X3) 
y = {—xi sin a + ui (x) cos a, —X2 sin a + U2 (x) cos a, U3 (x)) 



We show that 9Jt is also a "gradient" graph over x space. From Property 
2.2, we know that u {x'jX^) is a convex function in terms of x' for |x| < 1, 
or if necessary |x| < with depending only on m. From (j2.2p we also 
assume \D'u3 {x)\ = |(tti3,ii23) {x)\ < 1/2 for |x| < Vm- Then we have 



(3.2) 85 (m) \x — x*p > \x (x) — X (x*)|^ 

= (x' — X*') cos a + 



D'u{x',X3) - D'u{x',xl) 
+D'u{x',xl) - D'u {x*',xl) 



smo 



> 



(x' - X*') cos Q + { u (x^xg) - D'u{x*\xl) ] 
- I {D'u (x', X3) — D'u (x', X3)) sina|^ + la^s — 2:3 



+ F3 - 2:31 
2 n 



sma 

*|2 



> 



cos^ a I „/ 1 2 



X — X I + COS a sm a (^x — x ,Du[x, X3 j — D 



u[x*',xl)) 



■ sin^ a 2 L''n3 



>o 



<1 



> 



cos^ a 



(3.3) 



X — X I + (^1 — sm aj |X3 — X3I 

^ 1 I *|2 
> - X — X . 

- 4 I 



It follows that 9Jt is a special Lagrangian graph (x,Dn(x)) over a domain 
containing a ball of radius l/\/2 in x space. The Hessian of the potential 
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function u satisfies 



dy_ 

dx 



dy_ 

dx 



dx 
dx 



sma + nil cos a 
ui2 cos a 

cosQ + nil sin a 
ni2 sin a 




ni2 cos Q 
- sin a + n22 cos a 

Ui2 sina 
cos a + U22 sin a 




ni3 cos Q 
n23 cos a 
U33 

ui3 sin a 
U23 sin a 
1 



-1 



(3.4) 
and 

(3.5) det D'^u 



tan (^j — a) 



tan (^ — a) 



/IT 

tan — 

.4 



a 



m n_ 9 /vr 
-— 2 + tan(- 



,m— 11 



a]H: 



33 



,2ni-ll 



where the above abused notation l^r"*"^] also represents a matrix whose aU 
entries are analytic functions starting from order m— 1, and ()3.4p (13.50 follow 
from a simple calculation and the asymptotic behavior of D^u, (j2.2|) . We 
verify the following three properties for D^u. There exists a positive number 
'>^m,a depending only on m and a G [0, vr/4) such that for |x| < fm,a we have: 
Property 3.1. The determinant detD^n(x) is negative for small x 7^ 0, 
indeed 

det D^u (x) « — tan ^— — ; 
Property 3.2. The upper left 2x2 principle minor of the Hessian D^u, 



/IT 

2 tan - 

.4 



a) / > {D^u)' > 



tan I 



Property 3.3. The three eigenvalues Aj of the Hessian D^u satisfy 
9i = arctan Ai 
02 = arctan A2 
^3 = arctan A3 




tan(f-a) 



where " ~" means two quantities are comparable up to a multiple of con- 
stant depending only on m and a. Relying on (j3.5p . repeating the argu- 
ments for the estimate of A3 in Section 2, using (|3.3p and (j3.ip . we obtain 
Property 3.1. Property 3.2 follows from ([331). From ([331) dMI and the 
Lipschitz continuity of eigenvalues in terms of matrix entries, we derive the 
estimates for the first two eigenvalues in Property 3.3. In turn, noticing 

A3 = det D'^u/ ^AiA2^ , relying on both p.2p and ()3.3p we get two sided 

estimates of the last eigenvalue. 
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(3.6) 



> kI 



for X G Bp. 



Step 2. We proceed with the following proposition. 

Proposition 3.1. Let L = [x, Df) he a Lagrangian surface in = M.^ x\ 
with the smooth potential f over Bp cM.^, satisfying: 

Df (0) = 0, 

det D^f{x) <0 forx^Q, 

fu (x) fi2 (x) 

/2I (x) f22 (x) 

\D'f3{x)\ = \{f,^,f23){x)\<l, say 

Then £ can be re-represented as a graph {x,y) = ^x,Df (x)^ over open set 

n = Df [Bi^2p^ with X + v/^y = e~^^ [x + V^y) and f E (n) n 
C°° {n\ {0}). 

Proof of Proposition 3.1. Note that the U (3) rotation by 7r/2 is {x,y) = 
{y, —x) . This proposition really says that the map Df has a (unique) con- 
tinuous inverse $ = —Df. 

Step 2.1. We first prove Df is one-to-one on B^2p. Consider a coordinate 
change given hy t = (x) = (/i (x) , /2 (x) , X3). Then the Jacobian of ^ is 

fll fl2 /l3 

/21 /22 /23 (x) = det 
1 



(3.7) det Da;* (x) = det 



fll fl2 
/2I /22 



(x) > 0. 



Hence * is a local diffeomorphism on Bp. Note that is actually a distance 
expansion map. We have for all x, x* in Bp 

2 



* (x) - ^' (x*^ 



D'fix',xs)-D'f (x#',X3) 
+D'f {x*',xs)-D'f (x*',x* 



+ 



X3 -xj 



> 



i \D'f (x', xa) - D'f {x*', X3) I' - \D'f {x*\ X3) - D' f (x#', xf) 



+ 



-M, 

X3 — X3 



\ \D'f (x', X3) - D'f (x#', X3) -k{x' - X*') +k{x' - X*') 



X3 — X3 



> 



D'f {x*',x,)-D'f (x*', 
D'f (x, X3) - D'f (x*', X3) -k(x' - x#') , K (x - x#') 



>o 

,/fl|2 



+1|a.(x'-x#0| -2||i?7llioo(5,) 



X3 — X3 



+ 



X3 — X3 



k2 

> — 

- 2 



x' - x#' 



+ i-2||^7| 



L°°(Bp) 



X3 — X3 



(3.8) 



k2 
> — 

- 2 



X — X 



# 
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where we used ()3.6p . Thus \I' is a "global" diffeomorphism on Bp. 

We claim that the ^f-image of Bp, ^{Bp) D B^j^ . Otherwise, let be 

a boundary point of ^(-Bp) in BK^ . We know there exist a sequence of 

points Xi G i?p such that "^{xi) goes to and Xj goes to ^ Bp as i 
goes to infinity. If x* G -Bp, then ^ (x*) = by the continuity of But 
this is impossible because ^' (x*) is an interior point of ^'(-Bp) under the 
diffeomorphism of If x* G dBp, from (13. Sp . we have 

\t*\ = lim l^'(xi) - 0| > lim -^\xi - 0| = -^|x#| = -^p. 
i-s>oo i-j>oo -^2 V2 v2 

This contradicts G BK_ . 

Prom p.6p (only the upper bound), then 

|*(x) - ^'(x#)| < 1 1 DZ)'/ Hi- (Bp) \x - x*\ < K-^\x - x*\, 

if follows that ^'"^ is also a distance expansion map with a factor k. Apply 
the arguments above we get ^~^(B\ ) D Bio^ or B\ D "^fiBio). 

Now for the injectivity of Df on Bi^2p, it suffices to show that 

y{t) = Dfo^-\t) = {h,t2, /3(x(t)) 
is one-to-one in B\^^. Suppose that y (t) = y (t*) , then 



Note that 

dys (^1,^2,0 
St, 



tf = ti, q = t2, y3(t*) = ysW- 



1 
1 

^ ^ %3(tl,t2,g) 
9*3 



I(tl,t2,€) 



= detA {Dfo^-^) = det detA^ 

for {ti,t2,0 + 0, where we used and detA/(x) < for x / 0. 

It follows that the function yj,{t\,t2,^ is strictly decreasing in ^. Now 
Haiti, t2,tf) = y3{ti,t2,t3) implies tf = t^. This shows that y = Df o v]/-! 
is one-to-one. 

So far we have obtained the inverse function <I> = (Df)^^ onil. = Df (^Bi^2p 

Step 2.2. We prove Df is an open map from Bp to M^. Since the Jacobian 
det D^f (x) 7^ for x 7^ 0, Df is already a local diffeomorphism for x 7^ 0. 
It suffices to show that the image of an open neighborhood of in -Bp, under 
Df, contains an open neighborhood of in M^. Since ^' is a diffeomorphism, 
we only need to show this property for Df o vj/-!. Indeed we only need to 
consider the image of the ball B!^^ of radius 2r] centered at t = for a 
small r] > 0. According to Step 2.1, y{ti,t2,-) is strictly decreasing in the 
third variable. So 2/i„ = 2/3(0,0,7/) < and 2h+ = y3{0,0,-ri) > 0. By 
continuity of y = Df o \I'^^, there exists 77' G (0,r/) such that 2/3(^1, t2,^) < 



SINGULAR SOLUTIONS 



13 



/i_ < and y3{ti,t2, —rj) > /i+ > for |(ti,t2)| < f]' ■ Then by intermediate 
value theorem (for function 7/3 (ti,t2, •))' t^i^ open set 

{|(yi,y2)| < r/'} X < ya < C Z)/oi&-i({|(ti,t2)| < 1]'] X {ItsI < r/}) 
CDfom-\B%). 

Thus Z)/ is an open map. 

Step 2.3. Now = Df{B^2p) is an open neighborhood of y = 0, and ^ is 
continuous on by the openness of Df. Lastly we find a potential for the 
Lagrangian submanifold C now represented as (x, — ^> (x)) . Let 



[ -<^^{s)dsi-^^{s)ds2-<^^is)ds3. 
Jo 



Because {—^ (x)) is symmetric — {D'^f) ^ when x ^ and $ is bounded, 
7(5) is well-defined on Q. Further we know / G (Q) n C°° {0}) . 
The proof for Proposition 3.1 is complete. □ 



Remark. For the purpose of Theorem 1.1, we can replace (j3.6p by a 
weaker condition (j3.7p det [L*^/]' > 0. Consequently we have no estimate 
on the size of the existing neighborhood supporting the solution u in this 
section, then u"* for each single m and 0. The stronger assumption ()3.6p is 
designed for Theorem 1.2 where we need a uniform control with respect to 
e on the valid radius for the solutions u^, then . Lastly there is another 
argument for the openness of the particular map Du = Df, relying on the 
uniqueness of the pre-image of y = (which can be also derived from the 
distance expansion property at the origin ()4.3p ). instead of using the strict 
monotonicity property in Step 2.2. 

Step 3. Equipped with Property 3.1, 3.2, and ()3.4p . we apply Proposition 
3.1 to our function u (x) with x replaced by x, x replaced hy x, p = rm,a, 
and K = tan — a) /2. Then we get a new function u (x) , defined on 
an open neighborhood of x = 0. By Property 3.3, the three eigen-angles 
9i = arctan Aj of D'^u (x) away from the origin satisfy 



(3.9) 



?i=^i-f = -f-« + o(l) 
^2 = ^2-f = -f -a + o(l) 



~ \ + 0(\Dii(3:)r~^ 



2 tan(^-o) 

where the positive number 5m,a (x) is bounded from both below and above 
uniformly with respect to x, and 



-2a. 
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We verify u is still a viscosity solution to (jl.ip with = —2a across the 
origin. For any quadratic Q touching u at the origin from below, we have 
under the "diagonalized" coordinate system for D^u (0) 



D^Q < 



tan {—J — a) 

tan [—J — a) 

oo 



It follows that the eigenvalues A* of D^Q must satisfy 



TT 7T 71 

arctan AT < — — — a, arctan A, < — — — a, and arctan Aq < — . 

4: 2 



Then the quadratic satisfies 



arctan A* < —2a. 



i=l 



Observe that we can never arrange any quadratic touching u from above at 
the origin. Then there is nothing to check. When those testing quadratics 
touch the smooth u away from the origin, the verification according to the 
definition of viscosity solutions is straightforward. Thus u is a viscosity 
solution to (jl.ip with Q = —2a in a neighborhood of the origin. 

Step 4. Lastly we verify that the solution u is in fact (7i>i/(2ni--i) i^^^i j^ot 
C^'^ for any <5 > l/(2m- 1) in a neighborhood of the origin. The latter is 
easy. From Property 2.3 and ()3.3p . we see that 



(0,0,X3,I?n(0,0,X3)) = (0,0, X3, [x^] , [xi™] , -2m£52— i + [^2m]) 
It follows that 

1^3-01 |X3| 



\Du (0, 0, xs) - Du (0)1" (2me + [X3]) \xz 



S I ~ |(2m-l)<5 



00 



as X3 — for any 5 > 2m-i ■ This shows that u is not C^'^ for any 6 > 
l/(2m- 1). 
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Next we prove that u is C^''^/i'^''^ i) by the argument in [NV]. Observe 
that for i = 1, 2, 3 



\ui (x) - Ui {x*)\ 
|l/(2m-l) 



\X — X 



2m- 1 



Ui [X*] 



|=2m-l 



(5) 



< C (m) (2m — 1) sup luj ix] 



|l/(2m-l) 
i2m-2 



2m- 1 



:2m-l 



{x) 



= 2m-l 



\X — X^ 



Dui [x] I 



< C (m) sup \ Du (x] 

X 

< C (m) sup \ Du (x^ 

< C7 (m) , 



i2m-2 



i2m-2 



1 



\DiL [x) 



i2m-2 



where we used the fact that the scalar function ^^/(^m-i) ig (71/(2™- 1) ^j^i-j 
for the first inequahty, and (j3.9p for the third inequahty. 

Finahy by scahng u"* (x) = u (rx) /r^ with vahd radius r imphcitly de- 
pending on Tfi and. tli6 TfYi in Step 1 (We need to make this dependence 
exphcit and then to have a uniform control with respect to e on the valid 
radius for the solutions in Section 4. Our guaranteed valid radius goes to 
zero as m goes to infinity.), the desired solutions in Theorem 1.1 with each 
fixed O G (—^,0] are achieved. By symmetry, — are the sought solutions 
with phase 6 G [0, f ). 



4. Rotate to smooth solutions: proof of Theorem 1.2 

In this section, we create the desired family of solutions by another cor- 
responding families of U (3) rotations in on top those two in Section 3. 
For any fixed O G [0, |), let 47 = ^ — G > 0. We start the construction 
by taking small positive numbers e G (0,7) and solution u with fixed m in 
Section 1. 

Step 1. We take the U (3) rotation in Step 1 of Section 3 with a = ^ — ^. 
The valid radius of the rotation and the estimates of the Hessian D'^u are 
still valid. To prepare the final rotations in the last Step of this section, we 
require the following estimates of D'^u with eigenvalues Aj by shrinking the 
radius for x or |x| < re : 

' e\ = arctan Af = f- | + f + o(l)>7 
f4.n J 6'| = arctanA| = f -1 + ^ + 0(1) >7 
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where again the positive 6m,a (x) is bounded from both below and above 
uniformly with respect to x and e, further the above estimates and re are 
both uniform with respect to e. 

Step 2. Exactly as in Step 3 of Section 3, we apply Proposition 3.1 with 
p = re and k = tan (| — a) /2 to {x,Du) to get the potential u'^ with 

{i,DiL^) for S: e Du (^1^2^^) . It follows from i^J^ and ([SJ]) that 



(4.2) 



^i = -f-f + f + o(l)>7-f 
^2 = -f-f + f + o(l)>7-f 
~h = ^-\o(l)\ 



for |x| = \Du{x)\ < fe- We need to show this fe and the above o(l) 

terms are still uniform with respect to e, and also the |o(l)| term for ^3 
never vanishes when the input x does not vanish (actually this |o(l)| can 
be made explicit enough by (j4.3p l. All these can be seen from the following 
inequalities 

(4.3) S6{m)\x\ > \3:{x)\ = \Du{x)\ > 67 (m) \xf"'-'^ . 

Indeed we see the first inequality by recalling (j3.ip 



Du (x (x)) = y (x) = (cos a D'u (x) — sin a x' (x)) , 
(^D'u,us) (x'jXs) = Du{x) G C^, 



and (13. 2p . We have to work a little harder for the second inequality. Because 
of ()2.2p and a G (—§7) f — ^7) , the following convexity for function 



{x ) 



cos a u (X 



',303) 



sma 



is available 



cos a 



(a 



sm a 



cos a — sm a ^ 
> z / > 



for |(x',X3)| < re, where we shrink re if necessary. Then we get 

|y (x)|2 = I (D<3 (x') , na (x)) I' = I {Du',^ (x') - bx' + bx' , U3 (x)) | 
> \bx'\ + 2(^Du'^^ (x') — 5x',6x') + jus (x)|^ 



>o 



> 6^ Ix'P + Ins (x)l^ 
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where we set b = (cos a — sin a) /2 for simplicity of notation. In order to 
bound \x'\ from below, we use Property 2.3 to obtain 

\D'u {x',xs)f < 2\D'u {x',X3) - D'u{0,X3)f + 2\D'u{0,X3)f 



<C„|xf + |[ 



,2ml 1 2 



or 



2ml 1 2 



D'n(x)| - |[r2"] 



Hence 



\y{x)f > \Du{x)\'^ 



„2ml 1 2 



where we assumed the positive h'^/Cm < 1 without loss of generality. By 
virtue of Property 2.4, we get 

52 

|y (a;)|^ > — 

\2 







|^2m-lj2 _ 


[^2m] 1 



> 



(55 (m) 



|2m-l 



for |x| < re and small positive 5j {m,a) , where again we shrink tq if nec- 
essary. By (j3.2p we arrive at the second inequality of ()4.3p . 




Step 3. We make a final family of f/ (3) rotations in : i = e^^^^z. 
Again because U (3) rotation preserves the length and complex structure, 
m = (1,1)^2=) for |f| < fe still a smooth special Lagrangian submanifold 
with parameterization 

: COS e + Dvf (x) sin e 
■ sin e + Dn^ (x) cos e 

We show that OK is still a "gradient" graph over x space. From ()4.2p we 
know that the function (x) + ^ tan (f ~ 7) l^i] is convex. We then have 

X (x) — X (x*) = |(x — X*) cose + (^Dvf (x) — Dvf (x*)) sine] 

(x — X*) Jcose — tan — 7) sine] + 
+ [{D^L^ \x) - Du^ (i*)) + (5 - F) tan (f - 7)] sine 



> < 



+2 



I (x — X* ^ 
cose— 
tan (1 — 7) sine 



2 

(cos e — tan (f — 7) sin e) | + 



sine^x — x* 



(L»n^ (x) - Z)n^ {x*)) + 
(x - X*) tan (I - 7) 



= *|2 



X — X* COS e 1 



tane 



>o 

> — X — X 



tan 7 J ~ 4 

provided we take e G (0, 7) even smaller. It follows that the smooth 5Jt is a 
special Lagrangian graph ( x, Du'^ ( ^ ) ) o'^sr a domain containing a ball of 
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radius ^fe in x space. The eigenvalues Af of the Hessian D'^u'^ satisfy 
(4.4) 



arctan Af 



arctan A^ 
arctan A3 



l + f -e + o(l) 
l + f -e + o(l) 



|o(l) 



It fohows that vf' is smooth and satisfies 



arctan Af + arctan A| + arctan A' 



-G in Bi-- . 



Finally set 



{x) 



'¥4 



Observe that the gradients of the potential functions, or the heights of the 
special Lagrangian graphs are kept uniformly bounded with respect to e 
under the above three families of \J (3) rotations. Combined with ()4.4p . we 
obtain the desired family of smooth solutions to (jl.ip with n = 3 and fixed 
G G [0, 7r/2). By symmetry, —u^ are the other family of solutions to p.ip 
with n = 3 and fixed Q G (— 7r/2,0]. 



5. Minimal surface system: proof of Theorem 1.3 



In this section, we prove Theorem 1.3. Take the singular solutions iH^' 
from Theorem 1.1 with G = and m = 2, 3, 4, • • • . Let 

From Property 2.2 and Proposition 3.1, we see that 

\DU^ {y)\ = \D^u^ {y)\ 



1 



\Du^ {y)\ 



2m-2 ■ 



Here " ~ " means two quantities are equivalent up to a multiple of constant 
depending only on the dimension and m. Then we have 

1 



1 



Dm™{_Bi) \X 



(2m-2)p 



det [D'^u"' (y)] 



dy 



dx 



1 



\x 



(2m-2)p 



I 1 2m— 2 r J J 

\x\ dxidx2dx3. 



It follows that 



U'^eW^'PiBi) for any p < ^"^ ^ ^ but U"" W^'^ {Bi) . 

2m — 2 
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We next show that C/"* satisfies (II. 2|) in the integral sense, namely 

3 



for all <^ G j . This is because the integrand is everywhere 

except at the origin and we have the following bound on the integrand near 
the origin. Diagonalizing D'^u'^, we see that 



3 



3 

E 

i=l 



;i+A?)---(i+Ai)^5,ci>' 



<C{3,m) |Z)2n'"| \D<^\ =C{3,m) \DU"''\ \D<^\ G L\ 

where we used again the fact (]3.9p that two of the eigenvalues of D'^vT^ are 
bounded. The first part of the Theorem 1.3 is proved. 

The second part of Theorem 1.3 is straightforward if we take = Du^ 
with smooth solutions in Theorem 1.2 for any fixed © £ ( — f i f ) • 
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